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Abstract 

From the application point of view, it is important to have a good upper bound for 
the maximum norm of the inverse of a given matrix A. In this paper we will give 
two simple and practical upper bounds for the maximum norm of the inverse of a 
Nekrasov matrix. 
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1 Introduction 



Throughout this paper, we will use the following notations: 
C"(M") for the complex (real) n dimensional vector space, 
C"'"(M"'") for collection of sX\ n x n matrices with complex (real) entries, 
:= {1, 2, . . . , n} for the set of indices, 

n 

ri{A) = ^ |aij|, i E N for the deleted z-th row sum. 

jeiv\{i} 

Also, we define hi{A) recursively: 



Preprint submitted to Elsevier 



30 December 2011 



j=i nji j=i+i 



I ^ij I • 



We always denote hy A = D — L — U the standard splitting of a matrix A into 
its diagonal (D), strict lower {—L) and strict upper {—U) triangular parts. 

All classes of matrices, which are considered in this paper, are subclasses of 
if-matrices, more precisely nonsingular if-matrices. This class is well-know 
from many applications, and it can be defined by: 

Definition 1 A matrix A = [aij] G C^'^is called an H -matrix if its compari- 
son matrix [A) = [niij] defined by 



{A) = [m,,] e C"'^ m,, 



I (^ii I ) ^ J 



is an M -matrix, i.e., {A) ""^ > 0. 

A well-known property of iJ-matrices is given by the following theorem (see 

Theorem 1 Let A = [oij] G C"'" be a nonsingular H -matrix. Then 

\A-^\ < {A)-\ 

The most important subclass of iJ-matrices is the class of strictly diagonally 
dominant [SDD) matrices, defined by: 

Definition 2 A matrix A = [oij] G C"'" is called SDD matrix if, for each 
i E N, it holds that 

\aii\ > ri{A). 

Finally, let us recall the famous result of Varah [8] , which gives us an upper 
bound for SDD matrices: 

11^ ^lloo < r-r. . 7TTT- (1) 

min( - ri{A)) 

Obviously, this upper bound works only for SDD matrices, and even then it 
is not always good enough. As we see in above written formula, the smaller 
min(|ajj| — ri{A)) is, the worse estimation will be. For that reason it can be 



2 



useful to obtain new upper bounds, for a wider class of matrices which, as we 
will see, will sometimes work better in the SDD case, too. The class which 
will be chosen for that purpose is the class of, so called, Nekrasov matrices. 

Definition 3 A matrix A = [aij] G C"'",??, > 2 is called Nekrasov matrix if, 
for each i E N, it holds that 



It is known, see, for example [2], [7], that the class of Nekrasov matrices 
contains SDD class, while, on the other hand, it is a subclass of //-matrices. 



2 New estimations 



We will present two possibilities to estimate the infinity norm of the inverse 
matrix of a Nekrasov matrix. Although they look very similar, numerical ex- 
amples will show that each of them can be better than the other one. 



Theorem 2 Suppose that A = [aij] G C"'" is a Nekrasov matrix. Then, 

Z^{A) 



max 

1^ lloo ^ — 



-111 ^ \aii\ 

KiA) 



1 — max ■ 



ieJV \aii\ 
and 

rnax Zi{A) 
""^ ~ min(|an| - hi{A)) ' 

where 

Zi(A):=l, ^,(A):=|:M^,(A) + 1, zeN\{l}. (4) 

j=i Fiil 



In order to prove these estimations, we will start with the following lemma, 
proved by Robert in [B]. 

Lemma 1 Given any matrix A = [aij] G C"'", n > 2, with an ^ for all 
i G N, then 



h{A) = \a^^\ i\D\ - \L\)-'\U\e , 

L J I 

where e ^ is the vector with all components equal to 1. 



(5) 
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An immediate corollary of this Lemma is the following characterization of 
Nekrasov matrices, given by Szulc in [7]: 



Theorem 3 A matrix A = [ojj] G C"'", n > 2 is a Nekrasov matrix if and 
only if 

i\D\ - \L\)-'\U\e < e, (6) 

i.e. if and only if E — {\D\ — |L|)^-'^|^7| is an SDD matrix, where E is the 
identity matrix.. 

Proof of Theorem\^ First of all, notice that elements of matrix (|Z)| — |L|)^^|[/| 
are nonnegative, because matrix {\D\ — \L\) is an M- matrix. From inequality 
(|6]) we see that sum of all elements in each row is less then 1, therefore we can 
conclude that all diagonal elements are also less than 1. 

Assume that matrix ^ is a Nekrasov matrix. Then, from Theorem |3l we know 
that E — {\D\ — \L\)~^\U\ is an SDD matrix, so we can apply Varah bound for 
estimation of the infinity norm of its inverse matrix: 



\C II oo ^ 



max 



\cii\-ri{Cy 
where C := E - {\D\ - \L\Y^\U\. 

According to the fact that all diagonal entries of matrix {\D\ — \L\)^^\U\ are 
less than 1, we have: 



1 



m-\L\)-'\u\ 



^{c) = Y.[m-\L\)-^u\ 



and, therefore: 



n 

|cd-r.(C) = l-E[(l^|-|^l)-'l^l 



■n 

l-Y.[i\D\-\L\)-^\U\e 



K{A) 



Hence, 



\C lloo 



1 1 



1 ; — ;- 1 — max ■ 



, »eAf \aii\ 

Now, in order to estimate ||A^^||oo, it only remains to find a link between 
matrices and A~^. 

Since 

C={\D\-\L\)-\A), 
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it holds that 

and, see Theorem [H 



{A) = i\D\-\L\)C, 



p-ioo < ii(^)-ioo < iic-iooii(ii^i - mr'Woo. 

Finally, because \D\ — \L\ is an M- matrix, 

11(1^1 -|^l)-^lloo = 11(11^1 -|^|)-^elU, 

and, if we denote by y := {\D\ — \L\)^^e, then e = {\D\ — \L\)y, or, by compo- 
nents: 

i-l 

|aii|yi = 1, \aii\yi = 1 + '^\aij\yj, ieN\{l}. 

i=i 

Since \aii\yi = Zi{A), see (jl]), we get 



max 

ieN la,-. 



||(|D|-|L|)-i|U= Ibllc 
and the first upper bound (j2]) is proved. 

In order to prove the second upper bound, we will multiply matrix E — {\D\ — 
\L\)^^\U\ with diagonal matrix \D\ from the left side, and denote: 

B:=\D\-\D\{\D\-\L\)-'\U\. (7) 

Clearly, this multiplication does not change SDD property, so matrix B is also 
an SDD matrix, and we can, again, use Varah bound for estimation of the 
infinity norm of its inverse matrix: 



\B ^Wnc, < 



1 



max 



\bii\-ri{B)' 



Knowing that all diagonal entries of matrix (|Z}| — |L|) ^\U\ are less than 1, 
we have: 



ar 



{\D\-\L\)-'\U\ 



rm = j:K\[m-\L\r'\u\ 



. . 5 



and, therefore: 



\bii\ - niB) = \au\ - E K\ - \L\)-'\U\ 

3=1 



I ^ii I ^ii 



D\-\L\)-^\U\e\=\a,,\-K{A). 
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Hence, 

\\B 00 < max- — ■ , , ^, = — — ^ — ■ r-TTTT- 

i€N \aii\ - hi{A) min(|aii| - hi{A)) 

Now, in order to estimate ||y4^"'^||oo, it only remains to find a link between 
matrices and A~^. Since 

B=\D\{\D\-\L\)-'{A), 

it holds that 

(A) = {E- \L\\D\-^)B, 

and, see Theorem [H 

\\A-^\\oo < \\{A)-^\\oo < \\B-^\\oo\\{E - \L\\D\-\^\\^. 
Finally, in order to estimate IKE" — |L| IDI^"*^)^^ ||oo, we start with 



[E-\L\\D\-^)-^\\^ = \\{E-\L\\D\-^)-\\ 



00 ) 



which is true because ii^ — |L||Z}| is an M matrix. 

If we denote z{A) := {E - \L\\D\-^)-^e, then e = {E - \L\\D\~^)z{A), or, by 
components: 



zi(A) = l, z^(A) = l + J2p]z,{A), teN\{l}, 

j=l I'^jjl 

which is, precisely, the definition of Zi{A), given by (jlj). 
Therefore, 

\\{E - \L\\D\-')-'\\oo = \\z{A)\\^ = max Zi{A), 
and the second upper bound ([3]) is also proved. □ 

Both estimations, given by ([2]) and ([3]), can be used for bounding the infinity 
norm of the inverse of a Gudkov matrix, defined as a matrix for which there 
exists a permutation matrix P, such that PAP^ is a Nekrasov matrix. How- 
ever, before applying either ([2]) or ([3]), one can find the permutation matrix P, 
in order to ensure the first row of PAP^ to be strictly diagonally dominant. 



3 Numerical examples 



Keeping in mind that SDD matrices are subset of Nekrasov matrices, we will 
compare our two bounds with Varah's one, when it is applicable. For this 
purpose we will consider the following six matrices: 
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-7 1 -0.2 2 
7 88 2 -3 
2 0.5 13 -2 

0.5 3 1 6 



A, 



8 1 -0.2 3.3 
7 13 2 -3 
-1.3 6.7 13 -2 
0.5 3 1 6 



A. 



21 -9.1 -4.2 -2.1 
-0.7 9.1 -4.2 -2.1 
-0.7 -0.7 4.9 -2.1 
-0.7 -0.7 -0.7 2.8 



, ^4 



5 1 0.2 2 

1 21 1 -3 

2 0.5 6.4 -2 
0.5 -1 1 9 



A, 



6 -3 -2 
-1 11 -8 
-7 -3 10 



, ^6 



8 -0.5 -0.5 -0.5 
-9 16 -5 -5 
-6 -4 15 -3 
-4.9 -0.9 -0.9 6 



The upper bounds for their inverse matrices are given in the following table: 



Matrix 


class 


exact 1 1 A 1 1 oo 


Varah 


m 


m 


Ai 


SDD 


0.1921 


0.6667 


0.3805 


0.5263 


A2 


SDD 


0.2390 


1 


0.8848 


0.6885 


A, 


SDD 


0.8759 


1.4286 


1.8076 


0.9676 


A, 


SDD 


0.2707 


0.5556 


0.6200 


0.7937 


A, 


Nekrasov 


1.1519 




1.4909 


2.4848 


A, 


Nekrasov 


0.4474 




1.1557 


0.5702 



Obviously, each estimation ([2]) or (I3]) can work better than the other one. So, 
in general case, for Nekrasov matrices, one can take the smallest estimation 
of these two. 



For SDD case, it can be, also, useful to invest more time and calculations in 
obtaining these new estimations, since they can be significantly better (for 
example, see matrix A2). 
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